Elasto-capillarity at the nanoscale: on the coupling between elasticity 
and surface energy in soft solids 
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The capillary forces exerted by liquid drops and bubbles on a soft solid are directly measured using molecular dynamics simula- 
tions. The force on the solid by the liquid near the contact line is not oriented along the liquid vapor interface nor perpendicular 
to the solid surface, as usually assumed, but points towards the liquid. It is shown that the elastic deformations induced by 
this force can only be explained if, contrary to an incompressible liquid, the surface stress is different from the surface energy. 
Using thermodynamic variations we show that the the surface stress and the surface energy can both be determined accurately by 
measuring the deformation of a slender body plunged in a liquid. The results obtained from molecular dynamics fully confirm 
those recently obtained experimentally [Marchand et al. Phys. Rev. Lett. 108, 094301 (2012)] for an elastomeric wire. 



■1 Introduction 

As largely demonstrated in the last two decades, elasticity 
plays an important role in surface physics. Phenomena such 
as surface reconstructioni* 2 -, surface segregation—, surface ad- 
sorption 4 , elastic instabilities^, self assembly^,, and nanos- 
tructuration 8 of crystalline solids are directly induced by sur- 
face stresses. In parallel, and almost without any connec- 
tion, the elastic deformations of sheets and rods^"— as well 
as gels and elastomers induced by capillary forces have been 
evidenced and investigated—"—. It has remained unclear to 
what extent it is important to distinguish surface tension from 
surface stress for these elasto-capillary phenomena. 

The definition and properties of surface stresses can be de- 
rived from thermodynamics, atomistic models and mechanics. 
These approaches are complementary and should in principle 
be consistent with one another. The simplest situation is an in- 
terface between a condensed phase and its vapor. Let us con- 
sider an extensive quantity, the density of which varies across 
the interface over the scale of a few molecular sizes. At a 
macroscopic scale, the density can be seen as homogeneous 
on both sides of the interface. However, the extensive quantity 
then presents an interfacial excess. For instance, the free en- 
ergy presents an interfacial excess called the "surface energy", 
denoted y. Dividing now the total volume of the condensed 
phase and its vapor into two subsystems with a dividing plane 
normal to the interface, the mechanical force between the two 
subsystems also presents an excess quantity called the "sur- 
face stress", which throughout the paper we refer to as T. This 
is a force per unit length acting parallel to the interface, origi- 
nating from molecular interactions. 
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In the very particular situation where the condensed phase 
is an incompressible liquid, it can be shown from the virtual 
work principle that the surface stress and the surface energy 
are strictly equal, i.e. T = y. The surface stress and surface 
energy are then unified into a single name "surface tension", 
and it is common to address capillary problems using either 
the thermodynamic, or the mechanical route 34 ' 35 . In a solid, 
by contrast, the surface energy a priori depends on the strain 
in the bulk and yields an additional elastic contribution to the 
surface stress^ 3 ^. More precisely, the difference between the 
surface stress and the surface energy is the derivative of the 
surface energy with respect to the strain. This result is known 
as the Shuttleworth equation^ 3 ^, 



(1) 



where e is the bulk strain parallel to the interface. The sub- 
scripts refer to the phases i and j on both sides of the interface. 
Indeed, the Shuttleworth equation also applies when the inter- 
face separates two condensed phases composed of different 
molecules. This is highly relevant for wetting phenomena, for 
which one naturally deals with liquid-solid interfaces. Such 
interfaces present an excess free energy ygi that, according to 
(Q]), is different from the surface stress T^. Once again, the 
surface stress is the excess force parallel to the interface and 
can be measured at the edge of any control volume that in- 
cludes the interface. Importantly, the force Y$l is exerted on a 
subsystem composed two types of molecules, solid and liquid. 
From a molecular perspective, Y$l is the resultant of all types 
of molecular forces: solid-liquid, solid-solid and liquid-liquid 
interactions. 

The aim of this paper is to explore the difference between 
surface energy and surface stress in the case of soft solids 
that are partially wetted by a liquid. A paradigmatic example 
of this effect consists of a drop (or bubble) on a deformable 
soli d 17 i 18 i 22 - 24 i 28 ~ 33 , as shown in Fig. Q] This problem has 
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Fig. 1 Molecular dynamics simulations of a drop (a) and a bubble (b) on a deformable substrate, both with contact angle = 90°, hence 
Ysl = Ysv- ( a ) Drop on a soft substrate. Top: Snapshot, blue particles are liquid atoms, red particles are solid atoms. Bottom: Local 
displacement (red arrows) of the solid due to the presence of the liquid drop (shown by liquid isodensity contours), (b) Bubble on a soft 
substrate. Top: Snapshot, red and blue particles are the same as in (a), cyan particles are gas atoms. Bottom: Local displacement (red arrows) 
of the solid due to the presence of the bubble (shown by liquid isodensity contours). Note that the tangential displacement near the contact line 
is different between both situations: in the drop case the solid is pulled inwards whereas in the bubble case the solid is pulled outwards. 



recently been explored experimentally and theoretically, with 
contradicting interpretations emerging from microscopic and 
macroscopic descriptions 24 i 27 ~ 30 i 32 ] 33 ] 37 . Indeed, we expect 
a particularly strong manifestation of the difference between 
surface energy and surface stress near a three-phase contact 
line. On the one hand, the liquid-vapor interface is character- 
ized by a surface tension, i.e. Ylv = ^lv, and Young's law for 
the equilibrium contact angle 



Ylv cos 9 = ysv - Ysl 



(2) 



only involves the surface energies. On the other hand, the me- 
chanical equilibrium of the solid involves the surface stresses 
Ysv + Ysv and T SL ^ Ysl- 

In this work, we reveal the connection between elasto- 
capillary interactions at the nanoscale and the thermodynamic 
concepts of surface energy and surface stress. In Section|2] we 
perform Molecular Dynamics simulations of drops and bub- 
bles on soft substrates, for which we can accurately determine 
the elastic deformations and the liquid-on-solid forces. We 
find that the deformation below a drop is markedly different 
from the deformation below a bubble (Fig. [1}: we measure a 
tangential force oriented towards the liquid side, even for con- 
tact angles of = 90° for which Ysv = Ysl- Can one explain 



this tangential force by invoking surface stresses? In Section 
[3] we therefore develop a purely thermodynamic view, in the 
case of a plate partially immersed in a liquid, and compare 
this directly to Molecular Dynamics simulations. Finally, we 
conclude in Section 4 by relating the tangential forces to the 
difference between surface energies and stresses. 



2 Drops and bubbles 
2.1 Molecular Dynamics 

The aim of this paper is to investigate the basic mechanisms 
controlling elasto-capillary interactions in model situations. 
We have therefore performed Molecular Dynamics simula- 
tions using simple interactions for both the liquid and the soft 
solid, in a quasi-2D geometry. The simulations have been per- 
formed using the GROMACS software package^. The liquid 
consists of a Lennard-Jones fluid whose pair interaction po- 
tential is given by: 
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Here, ey and dij are the interaction strength and range between 
particle types i and j, respectively. The potential is cut off at 
5(1ll, where dn is the liquid atom size. The simulations are 
performed in the A^VT-ensemble (constant number of parti- 
cles N, constant volume V, and constant temperature T using 
the thermostat described in 3 ?). The solid consists of atoms 
placed on a cubic lattice of 25 layers, with harmonic springs 
connecting each atom to its neighbour and next-nearest neigh- 
bour. The lattice spacing a = O.Sdix, while all spring con- 
stants are taken equal with k = 38.5 • ei±/d\ L . There are no 
Lennard- Jones solid-solid interactions, hence the solid atoms 
only interact with each other through the harmonic springs. 
The solid (S) interacts with the liquid (L) through Lennard- 
Jones interactions. By varying the solid-liquid interaction, we 
explore drops of different equilibrium contact angles 40 . In the 
simulations for a gas bubble, we added gas atoms (G) that also 
interact according to a Lennard- Jones potential. The addition 
of gas atoms in the bubble case is required to prevent the bub- 
ble from collapsing immediately, as would be the case for a 
vapour bubble. The Lennard-Jones interaction parameters are 
given in Table Q] In the following paragraph we define the 
relevant dimensionless numbers. 

From a macroscopic view, the elasto-capillary deformations 
arise from a balance between surface tension y and the elastic 
modulus E (for simplicity of notation we use y for the liquid- 
vapor surface tension Jlv)- The ratio of these parameters y/E 
gives the elastocapillary length, which sets the scale of the 
elastic deformations. Our quasi-two-dimensional simulations 
have plane-strain conditions, in which case the relevant elas- 
tic modulus reads E = E/(l — V 2 ), where E is the Young's 
modulus and v the Poisson ratio. In terms of lattice parame- 
ters in our simulations, we find E = 15k/ (8a). To quantify the 
relative softness of the substrate, one can compare the elas- 
tocapillary length to the (liquid) atomic size djx, which gives 
the dimensionless quantity y/ (Edn). Whenever this quantity 
is small, the deformations are weak and one should recover the 
contact angles according to Young's law—. Here we measure 
the liquid-vapour surface tension 7 in a separate system using 
a Kirkwood-Buff integral over the stress-anisotropy near the 
interfaces y = J(pn ~ Pr{z))dz where, p^ is the (constant) 
thermodynamic pressure in the system, and pj the tangen- 
tial (relative to the interface) component of the stress-tensor 
which deviates from near the liquid- vapour interface^^ 2 -. 
We find, using the method outlined in ? to determine the lo- 
cal pressure that y = 0.78 en/djj^. In this work, therefore, 
the parameter 7/ (Edn) = 8.6 • 10~ 3 is indeed small, meaning 
that all elastic displacements are much smaller than atomic 
size (small elastic strains), and allows the use of linear elastic- 
ity theory. The typical length scale of thermal fluctuations in 
the solid, y/kgT/k, compared to the elastocapillary deforma- 
tion described before provides another dimensionless quantity 

g^=y = 0.18. The smallness of this parameter shows that 



Table 1 Lennard-Jones Interaction parameters for the 
MD-simulations for Liquid, Solid and Gas atoms. With these valus, 
the liquid-vapour surface tension is 7= 3.1 • 10~ 2 J/m 2 . 
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0.34 


SL 


varied 


0.34 
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GG 
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0.5 


SG 


0.004 


0.34 


LG 
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0.42 



the thermal fluctuations in the system are much larger than the 
deformation due to capillary forces. Still, as we will show, 
the displacement field can be measured very accurately in the 
simulations after averaging over time. Finally, we note that 
in Si-units, the chosen material properties correspond to real 
materials at T = 300K, E = 11 GPa, and 7 = 3.1 • 10~ 2 J/m 2 . 
Surface tension coefficients of simple liquids typically lie be- 
tween 0.02 J/m 2 (ethanol) and 0.07 J/m 2 (water). Young's 
modulus of crystalline solids is typically around E = 100 GPa 
while it can be much lower for elastomers (between 10^2 and 
10 1 GPa for rubber) and gels (down to 1 kPa) whose elastic- 
ity is entropic. 

The small strains in the solid (smaller than thermal fluctua- 
tions) are measured by calculating the time-averaged displace- 
ments (relative to the center of mass of the droplet or bubble) 
of the solid atoms compared to a base state, obtained from a 
simulation of the same solid in vacuum. As we are interested 
in the influence of the liquid on the solid, this procedure al- 
lows us to exclude effects associated with the presence of a 
solid- vacuum interface. The contact angle of the droplet and 
of the bubble are measured by determining the time-averaged 
density field of the liquid and the position of the Gibbs inter- 
face. We refer the reader to our previous work 41 ' for technical 
details. The circular fit to this liquid-vapor/gas interface is ex- 
trapolated to the solid, which provides the contact angle. 

2.2 Elasto-capillary deformations 

To illustrate that surface energy is not sufficient to characterize 
elasto-capillary deformations, we first consider a case where 
Ysv = ysL, sucn that the contact angle of the liquid is close 
to 90°. The contact angle can be adjusted by tuning only the 
Lennard-Jones interaction e^u while keeping the liquid pa- 
rameters (ei±, <7fx, <Jsl) fixed. This way, the liquid properties 
are unchanged except for the interaction with the solid. Note 
that, in general, the surface stresses will be different from the 
surface energies, and thus Ysv 7^ Tsl- 

Figure Q] shows the elastic deformation in the solid below 
a liquid drop (panel a), and below a bubble filled with gas 
(panel b). Since 9 = 90°, the shape of the liquid-vapor in- 
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terface, characterized by the iso-density profiles, is very sim- 
ilar in both cases. By contrast, the elastic deformations are 
markedly different, as can be seen from the vector field (red 
arrows): while below the drop one observes a tangential dis- 
placement towards the center of the drop, the displacements 
below the bubble are oriented outwards. This surprising out- 
come has important consequences. The drop and the bubble 
are perfectly symmetric from the point of view of the surface 
energies, since jsv = Ysl in this case. Yet, this symmetry is 
not reflected in the surface displacements: the deformations 
are not invariant under an inversion of the phases. Instead, the 
solid is always pulled towards the liquid side of the contact 
line (not only for the case = 90°). Therefore, an elasto- 
capillary description based on constant surface energies (i.e. 
on surface tensions) is not sufficient to describe the elastic de- 
formations below a drop or below a bubble. 

2.3 Capillary traction and contact line force 

Before turning to a fully thermodynamic description of elas to- 
capillary deformations in Sec. [3] we first quantify the capillary 
liquid on solid forces at the nanoscale. Note that for a com- 
plete description of the bulk deformations one requires both 
liquid on solid forces and any solid on solid forces that are 
present in the surface layer. However, since the asymmetry is 
caused by the liquid on solid forces we first quantify the liquid 
on solid forces only. The common feature of the deformation 
below the drop and the bubble is that, below the contact line, 
the tangential deformations are oriented towards the side of 
the liquid phase (Fig. [TJ. This is consistent with the predic- 
tions of the Density Functional Theory in the sharp interface 
approximatio n 29 ! 35 , which is based on a microscopic descrip- 
tion of the interactions. The mechanism for this asymmetry is 
that the long-ranged attraction by the liquid molecules creates 
a resultant force on the solid that is biased towards the liquid: 
the solid is more strongly attracted by the phase of highest 
density. This resultant force ultimately determines the elastic 
deformations, and is responsible for breaking the symmetry 
between drops and bubbles in Fig. Q] 

In Molecular Dynamics, we can of course directly quantify 
this effect by measuring the time-averaged forces that all liq- 
uid molecules exert on the solid molecules. Since we wish 
to reveal the capillary traction (force per area) that the liquid 
induces on the solid, we integrate over the vertical depth of 
the solid. In practice, the direct liquid-solid interaction only 
acts in the superficial layers of the solid, due to the short range 
of Lennard-Jones (van der Waals) interactions (~ r~ 6 ). The 
resulting capillary traction below a drop is shown in Fig. |2] 
First, the curve in Fig. [2^ shows the normal traction, a,,, of 
the liquid on the solid in the case = 82°. As expected, we 
observe a large traction that is localized near the two contact 
lines: this corresponds to the "pulling" action of the contact 



line. The width of the peak is a few molecular sizes and re- 
flects the width of the liquid- vapor interface. At the center of 
the drop one observes a slightly negative traction, correspond- 
ing to the Laplace pressure in the bulk of the drop. As the drop 
is in equilibrium, this Laplace pressure perfectly balances the 
upward stress at the contact line: the total liquid-on-solid nor- 
mal force is zero. Using macroscopic thermodynamics we can 
estimate the normal force exerted by the liquid on the solid in 
the vicinity of the contact line. For drop sizes that are much 
larger than the width of the peak in a„ (which coincides with 
the thickness of the liquid-vapour interface), we can separate 
the capillarity forces into a contribution per unit area, P = y/R, 
and a perfectly localized force per unit contact line /„. Using 
that the width of the drop 2R sin 6 and that the total force van- 
ishes, one predicts that the strength of the force on the solid 
near the contact line should be /„ = ysin (per unit contact 
line). We can test this macroscopic prediction in our simu- 
lations. First, we determine the liquid-vapor surface tension 
y from an independent calibration, as described in Sec. 12.11 
Then, we integrate the normal stress over the peak located 
around the contact line, yielding the total normal force per 
unit contact line /„. The result for different contact angles 6 is 
shown in Fig . |2j: . The results of MD simulations are consistent 
with a contact line force in the normal direction /„ = ysinf?, 
shown in solid line. 

Similarly, we can determine the tangential capillary trac- 
tion on the solid, denoted a t . The result is shown in Fig. ^p. 
We find a positive traction at the contact line located on the 
left (i.e. pointing towards the right), and a negative traction 
at the contact line located on the right (i.e. pointing towards 
the left). Indeed, we identify an effective liquid-on-solid force 
that is oriented towards the interior of the drop, i.e. into the 
liquid phase. This tangential force is the reason why the elas- 
tic displacements point towards the interior of the drop. Again, 
we quantify the total force exerted on the solid near the con- 
tact line from the integral of the peaks, The resulting tan- 
gential force per unit contact line is shown in Fig. |2jl, as a 
function of the contact angle. We observe a nonzero inward 
tangential force for all angles. This can be understood from 
the left-right symmetry breaking below the contact line: the 
solid atoms are attracted much more strongly by the high- 
density liquid phase than by the low-density vapor phase. As 
the strength of the solid-liquid interaction is directly quanti- 
fied by the work of adhesion, y+ jsv — Ysl, one expects below 
a liquid/vapour/solid triple line^o 

ft =y+75V-75L = y(l+cos0) . (4) 

This equation is shown as the solid line in Fig. [2jl: it indeed 
captures the features of the tangential force, which is always 



* In case the vapor is replaced by a second, immiscible liquid, there will be sec- 
ond "work of adhesion" contribution. This yields another tangential force on 
the solid that is biased toward the second liquid, and weakens the asymmetry. 
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Fig. 2 Capillary traction: the liquid on solid forces per unit area of the substrate measured in molecular dynamics of a droplet on a substrate 
(9 = 82°). (a) Normal component of the force per area by the liquid on the solid, rr„. The two peaks correspond to the contact lines, where the 
solid gets pulled up. The region between the peaks corresponds to the liquid-solid interface, where the solid gets pushed down due to the 
Laplace pressure, (b) Tangential component of the force per area by the liquid on the solid, a,. There is only a force near the contact lines, and 
the force is directed towards the interior of the droplet, (c ) Normal component of the total force by the liquid on the solid due to the contact 
line at varying 9. The solid line corresponds to sin 9. (d) Tangential component of the total force by the liquid on the solid due to the contact 
line at varying 9. The solid line corresponds to 1 +cos 9. The low values for small contact angles 9 < 40° are likely due to finite size effects 
in the simulations. 



positive, i.e. oriented towards the liquid side, and nicely de- 
scribes the magnitude and trend with the contact angle. This 
also explains the difference between the deformations below a 
drop and a bubble. 

In conclusion, our simulations clearly demonstrate the ex- 
istence of a tangential capillary force exerted on the region 
of the solid below the contact line. This force has a strong 
influence on the elastic deformation below a drop or bub- 
ble. This effect is usually ignored in the literature on elasto- 
capillarity 17 i 18 i 22 i 28 i 30 i 32 , likely due to the fact that in the case 
of an incompressible liquid (or solid) this tangential force is 
exactly balanced by solid on solid forces and therefore is not 
transmitted to the bulk. However, in general for v ^ 1/2, 
this is not the case and the tangential force (which is always 



pointed towards the liquid phase, see Fig.[2ji) needs to be taken 
into account when considering the elastic deformation. 

3 Partially immersed solid 

From the preceding section is it clear that the elastic deforma- 
tion below a contact line results from the detailed interactions 
(capillary and elastic) at the nanoscale. We will now address 
the problem in a macroscopic framework, where we relate the 
elastic displacements to purely thermodynamic concepts. In 
particular, the goal is to express the microscopic interactions 
discussed in Sec. 2 directly in terms of the excess quantities 
Yu and Ty. For this, we consider an long elastic plate that 
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is partially immersed in a liquid - see Fig. [3] This geome- 
try was studied before experimentally using a thin elastomeric 
wireM. The experiment revealed that the elastic strain in the 
"wet" part of the solid was very different from the strain in 
the "dry" part. Here we analyze this geometry using the ther- 
modynamic concepts of surface stress and then compare it di- 
rectly to Molecular Dynamics. 

The central result of this section is that the vertical strain 
above the contact line, e+, couples to the surface energies, 
while the strain below the contact line, e_, is determined by 
the surface stresses: 

WE - 

— e+ = Ysv-Ysl (5) 
WE 

— £_ = T SV -T SL . (6) 

Here W is the width of a two-dimensional elastic plate, E is 
the elastic modulus, while the reference state for the strain is 
the completely dry solid. By determining the strain inside the 
plate, one thus directly measures the difference between sur- 
face energies and surface stresses. This method will be applied 
in our Molecular Dynamics simulations. 



To derive (O and @, we consider the free energy of an elas- 
tic plate that is partially submerged in a liquid bath with the 
contact line at z — (Fig. [3]). We will assume that the elasto- 
capillary length y/E is much smaller than the plate width W, 
which in turn is much smaller than the total plate length L: 

I < W < L . (7) 
E 

The first assumption allows one to consider the solid-liquid 
and solid-vapor interfaces to be flat with respect to other scales 
in the problem, as can be seen from the zoom near the con- 
tact line in Fig. [3J?. The second assumption implies that, apart 
from the region directly below the contact line of width W, the 
strain is homogeneous and characterized by a constant value 
for e = du z /dz, where u z is the vertical displacement field. We 
denote the strains above and below the contact line as e + and 
e respectively; the contribution to the energy of the region di- 
rectly below the contact line is sub-dominant by a factor W jh. 
Outside the contact line region the strains can be considered 
small, since e± ~ y/EW <C 1. Under these assumptions, we 
arrive at the following free energy functional of the system 
(per unit length): 



3.1 Thermodynamics: strain above and below meniscus 
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Fig. 3 Partially immersed plate of width W, held at equilibrium by 
an external force F ext . (a) The plate is partially wetted, with the 
contact line located at z = 0. There is a homogeneous strain above 
and below the contact line: e+ and £_, respectively, (b) Zoom 
around the contact line. Near the contact line exists a transition 
region of the strain from e + to e_ . The transition occurs over a 
length ~ W. 



& = 2 



dz 



y(l+h' 2 ) l/2 + y SL 



2 / dzysv 

-■cl 



d . 1 



+W I dz-Ee z _ 



W 



1 



dz -Eei 



(8) 



Here a and b denoted the bottom and top positions of the plate, 
but note that the elastic energy should be taken over a domain 
of fixed length Lq =bo~ fln> with reference positions «o an d 
£> ii. We allow for variations of the shape of the liquid-vapor 
interface h(z), the position of the contact line z c u the top and 
bottom positions of the plate b and a, and the elastic strains 
E± . The equilibrium conditions follow from 



S&=F ext Sb, 



(9) 



which equates the change in energy to the work done by the 
external force. Note that the positions b and a are linked by 
the constraint 



b — a = Lq ■ 



dz £■ 



(10) 



First, we derive the equilibrium conditions for the liquid, 
by considering variations of the liquid-vapor interface 8h(z), 
with 8b = 8e+ = 8e = 0. From geometry near the contact 
line, this implies a variation of the contact line position ac- 
cording to 8h(z c i) = —h'(zci)8z c i- One thus obtains 
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The integral expresses the Laplace pressure condition for the 
liquid-vapor interface. The terms ~ e 2 arise from the fact that 
a variation of z c i does not affect the contact line zone, but just 
gives an exchange of the elastic energies of the dry and wet 
parts respectively (similar to the exchange of surface energies 
Ysl — Ysv)- Collecting the terms from the boundary condition, 
using 8h(z c i) = —h'(z c i)8z c i, one finds the condition for the 
equilibrium contact angle 



also anticipate that the terms of order EWe± can be neglected 
in the following steps. 

The variation of the plate energy must be balanced with the 
work done by the external force 

S&p = F ext 8b = F ext 8bo + F ex ,b 8e + . (16) 

Combining (JT3J and ( fTol l. gives the three equilibrium condi- 
tions. The variation of bo gives the familiar expression for the 
external force needed to hold the plate, i.e. 

Fa =2(Ysv-Ysl)- (17) 
Using this expression, we find that the strains must follow 

-dYsv 



EWe+ 
EWe- 



de 



2y< 



SL 



(18) 
(19) 



cose = + 

y 2 y ' 



Ysv - Ysl 
Y 



\EWJ 



(12) 



where we replaced cos = 1/(1 + h' 2 ) 1 ^ 2 . This shows that for 
Y/EW <C 1 one recovers Young's law for the liquid contact 
angle with respect to the undeformed solid. 

Next, we explore the elastic degrees of freedom of the plate. 
For convenience, we now choose the contact line position as 
the reference altitude: z c i = 0. Using Eq. (8), one writes the 
free energy & p of the plate and its interfaces with the liquid 
and the vapor: 

& p = 2bYsv (e+ ) - 2a 7 sL (£- ) + \ WE (b Q £ 2 + - a Q e 2 _ ) . (13) 

Here we made explicit the strain dependence of surface ener- 
gies, which is necessary for solid interfaces. Due to the rela- 
tions 



The reference state for the strain, £o, will be the plate in con- 
tact with the vapour only. This reference can be derived in a 
similar way by minimizing the free energy of a plate in vapour 
(i.e. replacing Ysl by Ysv i n either of the above equations): 



EWeo = -2^-2Ysv 



(20) 



The final step is to subtract the reference state from (TT81 and 



( fT9l ). and substitute the Shuttleworth equation Y ;; - = + Yij 



£+ — £+ — £q 



£ - £o = 



EW 



(r w -T rf ) 



Indeed, this is the result anticipated in (0 and ©. In conclu- 
sion, the immersion of a slender body allows one to determine 
accurately both surface energies and surface stresses—. 



b = b (l+e+), fl = ao(l+e_), a = b -L , (14) 3.2 Molecular Dynamics 



there are only three independent variables. We choose here to 
parametrize the problem using bo, £ + and e . From this we 
can write the total variation: 



1 



8^ 



8bo 
+b 8e 
—aoSs 



Ysv - Ysl + 



Ysv 



Ysl 



1 

2 

dYsv 



de 

dYSL 

de 



EW (4 - 

EW 
f 

2 

EW 



(15) 



where, after the variation, we replaced b — bo and a = «o ow- 
ing to the smallness of the strains. For the same reason, we 



We now test this method in a MD simulation of an immersed 
plate. Figure |4] shows a snapshot of the simulation, where the 
solid plate is partially immersed into a liquid reservoir. The 
material of the plate is the same as that used in the previous 
drop and bubble simulation, except that it is stiffer: k p = 10k, 
such that £p = 1 10 GPa whereas the liquid surface tension is 
unchanged: y= 3.1 ■ 10~ 2 J/m 2 . The container that holds the 
liquid consists of the same material as the plate, except that the 
container is not allowed to deform by fixing the atoms to their 
initial positions. To avoid a curved liquid meniscus, and hence 
a difference in pressure on the wetted and dry regions of the 
plate, we consider again 9 = 90°. In this case Ysv = Ysl an d 
according to (0 we expect e + = 0. Before each simulation 
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Fig. 4 Snapshot from an MD-simulation of a plate (W /a = 26 
atoms wide) partially submerged in a liquid bath. The liquid-solid 
interaction energy cj/, was chosen such that 6 = 90° . Comparing 
the time-averaged z-position of the solid atoms to the time-averaged 
position of the reference system (plate in vacuum), allows for the 
local displacement ;( ; (z) to be measured. 



we equilibrated the solid in a vacuum to have a well-defined 
reference state. 

Figure shows an example of the vertical elastic displace- 
ment, u z (z), that is induced after immersion of the plate. The 
slope of this curve directly gives the strain e = du-jdz. In- 
deed, we observe very different strains above and below the 
contact line, which allows for a determination of e + and e_ . 
In this example the top part of the plate is hardly deformed, 
as expected from (0 for this situation where Ysv = Ysl- By 
contrast, the lower part of the plate displays a negative strain, 
e_ < 0, corresponding to a compression of the bottom part 
of the wire. Using ©, this reveals a difference in surface 
stresses, Ysv 7^ ^sl, despite the equality of surface energies. 

To quantify the difference Ysv ~ ^"sl, we repeated the simu- 
lations plates of various widths W. The resulting £ + and e_ are 
plotted as a function of 1 /W (FigO. Above the contact line 
we indeed find a vanishing strain e + for all plate thickness, 
within the error bars of the simulation. Below the contact line 
we observe a nearly linear dependence on 1 /W, as predicted 
by (|6). This confirms that our system is large enough to apply 
thermodynamics and continuum elasticity. The slope of the 
curve is measured at = — (1.4±0.3) • 10~ 3 di±. This 

slope can be related to the surface stress via Eq. (|6} and we 
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Fig. 5 MD measurement of the relative displacement u~ (z) for pa 
late width W = S.Sdn. The derivative of this slope gives the strain 
e(z). The contact line is located at z = 0, where a clear jump in e(z) 
is observed (from e_ to £+, see also Figs.[5]and|4]) This signifies an 
imbalance in surface stresses (Y^ ^ Yjy) even though Ysl = Ysv f° r 

e = 90°. 

find, for this specific solid and liquid pair: 

"Csv - r SL - (Ysv - Ysl) = Vsv ~T SL = - (0.8 1 ± 0. 17) y . 

(21) 

This surface stress difference is clearly not a negligible effect: 
it is of the same order as the liquid-vapor surface tension. 

4 Conclusions 

We have shown from thermodynamic considerations that the 
elastic deformation of a partially wetted solid crucially de- 
pends on the difference between surface stress and surface 
energy. This result is confirmed using Molecular Dynamics 
simulations, revealing how deformations emerge from inter- 
actions at the nanoscale. To complete the picture of elasto- 
capillary interactions, we finally give a purely mechanical in- 
terpretation of our findings. 

Describing the partially immersed wire of Fig. fusing con- 
tinuum elasticity, a discontinuity of strain implies a discon- 
tinuity of stress across the contact line, in the direction par- 
allel to the solid interface. This means that the contact line 
region must exert a tangential force ff on the bulk elastic ma- 
terial—. The magnitude of the tangential force experienced by 
bulk elasticity is proportional to e + — £_, and therefore reads 

ff = (Vsl ~ Vsv) - (Ysl -Ysv)- (22) 

This residual force accounts for all interactions that are trans- 
mitted across the surface layers to the bulk elastic, including 
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Fig. 6 MD measurements of the vertical strains £ + (circles) and e_ 
(squares) as functions of the inverse plate width dn/W . Within 
error, there exists no strain above the contact line (£+). This is 
expected for = 90°. Below the contact line, however, the solid is 
compressed due to an imbalance of the surfaces stresses at the 
contact line (Yj^ 7^ Y^y). The slope of this curve quantifies the 
magnitude of this imbalance, Eq. l(6j. 
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Fig. 7 Stresses acting on the circular control volume around the 
contact line. Note that for solids generally T^x 7^ 7sx> hence there 
exists an imbalance of the interfacial stresses in both the normal and 
tangential directions. This imbalance is counteracted by elastic 
stresses in the solid, along the thick gray dashed line. 



the solid-solid interactions. It is therefore important to dis- 
tinguish ff from f t measured in Fig. [2] while the latter only 
included the liquid-on-solid forces, the surface stress captures 
the total excess surface force and includes all superficial inter- 
actions. 

Figure [7j shows how the residual force ff arises due to 
the imbalance of surface stresses in the vicinity of the con- 
tact line. When discussing the forces near the contact line, it 
is absolutely critical to explicitly specify the material system 



to which the forces are applied: a different choice of control 
volume will lead to different force s 35 ' 43 . Here we consider 
a macroscopic control volume that includes the three-phase 
contact line, as indicated by the dotted circle. As this includes 
the three interfaces, one can directly represents the surface 
stresses Y$y, Y$i and T = y as indicated by the solid arrows. 
Interestingly, the equilibrium contact angle does not involve 
the surface stresses of the solid, but rather the surface ener- 
gies ysi and y$v- This is not inconsistent with Fig. [7] because 
Young's law represents an equilibrium (a minimal free energy 
from the thermodynamic perspective and a balance of forces 
from the mechanical point of view) inside the liquid only, and 
thus requires a different control volume that does not include 
the solid — . As a consequence, the surface stresses in Fig. |7] 
do not balance in the direction parallel to the solid, but yields 
a nonzero tangential force ff = ycos 9 +Tsl — Xsy, in agree- 
ment with (f22b . Similarly, the surface stresses yield a resultant 
normal force 

ff = ysm9. (23) 

To restore mechanical equilibrium inside the control volume, 
both ff and ff must be balanced by elastic stresses that are 
exerted along the along the grey dashed circular section in 
Fig- El The tangential component vanishes only when the sur- 
face energies and surface stresses are equal. 

A similar observation can be made for very soft solids, for 
which the solid deforms into a "cusp" shape with a solid angle 
Q s < 7T 28 i 3Q ' 37 . The cusp develops when y/Ed LL > \2L, and 
thus corresponds to cases much softer than in our Molecular 
Dynamics. Figure [8] shows the surface stresses near the con- 
tact line on such a strongly deformed solid. For given value 
of 0s, the equilibration of the liquid angle 6l involves only 
the surface energies y$L, ysv an d n °t the stresses Ysl^sv- I n 
general, the liquid equilibrium will therefore not coincide with 
the balance of surface stresses on the circular control volume 
in Fig.[8p 7 : residual elastic stress will arise whenever surface 
energies differ from surface stresses. 

The key parameter for wetting of soft materials is thus the 
difference between surface energies and surface stresses, as in 
equation (l22l . Experimentally, the strain discontinuity across 
the contact line for the partially immersed wire gives direct 
access to this difference. The recently suggested method to 
determine the surface stress from contact angles assumes a 
perfect balance of surface stresses—, and therefore incorrectly 
assumes that elastic stress can be ignored for the balance that 
determines the contact angle. This method, therefore, ap- 
plies only when there is no difference between surface ener- 
gies and stresses. Theoretically, we can now put upper and 
lower bounds on the tangential force. Clearly, the simulations 
in Fig. [TJ show that the tangential force originates from the 
breaking of left/right symmetry near the contact line, biased 
towards the side of the high-density liquid. The maximum 
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Fig. 8 Stresses acting on the circular control volume around the 
contact line, in case of strong solid deformation (y/Ea 1). The 
liquid equilibrates at an angle 81 that is a function of the solid angle 
Bs and the surface energies Ysl, Ysv- Since for solids generally 
Y$x Ysx> the liquid equilibration does not coincide with a balance 
of the interfacial stresses. This imbalance is counteracted by elastic 
stresses in the solid, along the thick gray dashed line. 

possible residual force ff should therefore be the liquid-on- 
solid force 7+ Ysv — Ysl- This maximum arises whenever the 
solid-solid interactions in the surface layer do not counteract 
this effect and the full tangential force is transmitted to the 
bulk substrate. This was referred to as the "vectorial force 
transmission model" in previous work—. Another extreme 
limit corresponds to the "normal force transmission model", 
for which the surface layers completely screen out any tan- 
gential stress. This is the case, for example, when the sub- 
strate is another liquid (i.e. an oil drop floating on liquid). A 
liquid can of course not sustain any shear, which means that 
the symmetry -breaking of interactions is counteracted by self- 
interactions inside the liquid substrate. This once more agrees 
with (|22| |. since for liquid-liquid interfaces T, 7 - = "fu. Hence, 
we conclude 

<ff < Y+Ysv-Ysl, (24) 

or equivalently 

Ysl - Ysv < ^sl - ^sv < Y- (25) 

Future work should be dedicated to determining how the dif- 
ference between T and y exactly depend on the material di- 
mensionless parameters: the solid Poission ratio the ratio of 
the elasto-capillary length to the atomic size and the ratio 
of the elasto-capillary length to the thermal length. The re- 
sults from our Molecular Dynamics are very close to the up- 
per bound: the surface stress difference was found slightly 
smaller than the liquid-vapor surface tension 7 A similar con- 
clusion can be drawn from the experimental results reported 
in—, where an elastomeric wire that was partially immersed 
in a liquid. The symbol T used in this previous work can now 



be identified with Y$l — Tsy, which was found identical to 
Y within experimental uncertainty—. To further explore the 
difference between surface energy and surface stress, it would 
be interesting to directly measure the tangential displacements 
inside soft substrates, e.g. using confocal microscopy—, and 
compare the deformation below a drop and a bubble. 
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